Identification of Variable Stars in CoRoT's first main observing field (LRcl) 
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ABSTRACT 

The COROT space mission will monitor several target fields for up to 150 days to 
perform asteroseismology and to search for extrasolar planets by photometric transits. 
Variable stars in the target fields are important objects for additional scientific studies 
but can also disturb the search for planetary transits. A variability characterization of 
the target fields prior to COi?OT observations is therefore important for two reasons: to 
find interesting variable stars to monitor further and to make an analysis of the impact 
of the variable stars on detecting extrasolar planet transits with COROT. 

The Berlin Exoplanet Search Telescope (BEST) is a small wide-angle telescope ded- 
icated to high-precision photometry. It has observed a 9 deg^ field of view centered at 
(a,5)=(19''00™00!0,-F00°01'55'.'2) (J2000.0) over 98 nights to search for variable stars in 
the surroundings of the first long-run target field (LRcl) of the COROT space mission. 
In this data set we identified 92 periodic variable stars, 86 of which are new discover- 
ies and 6 of which are known from the General Catalogue of Variable Stars (GOVS). 
For five of the GOVS stars, variability could not be confirmed. Forty-three of the 92 
detected periodic variable stars are identified as eclipsing binaries. We have evaluated 
the completeness of our survey for eclipsing binaries by comparing it to the expected 
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fraction of eclipsing binaries based on i^/PP^di^CO^ observations. From this evaluation 
we show that the BEST data set presented here has a completeness of 20% - 30% for 
periods longer than 1 day and is complete relative to HIPPARCOS for short-period 
binaries. 

Subject headings: methods: data analysis - binaries: eclipsing - stars: variables: gen- 
eral. 



1. Introduction 



Small wide-angle telescopes such as the Berlin Exoplanet Search Telescope (BEST: lRauer et al 



2004 ) offer a unique possibility of identifying stellar variability in a large number of stars. Apart 
from giving us important information on the individual targets, these kinds of wide-field data also 
contain information on the statistical distribution of variable stars, which can be used not only 
to understand the structure and evolution of the G alaxy, but a lso to analyze the expected false- 
alarm rates for transiting extrasolar planet searches (jBrownll2003l ) . For space missions searching for 
extrasolar planets by photometric transits, knowledge of the variable stars in the target fields helps 
reduce the number of false alarms, thereby reducing the amount of follow-up observation needed. 
In addition, the variable stars provide interesting targets for additional science. 



The COP OT space mission (IBaglin et al 



2002 



Borde et al.ll2003l ) , launched on 2006 December 



27, will observe a number of target fields for up to 150 days in order to perform asteroseismic studies 
of the target stars and to search for extrasolar planets by photometric transits. The observation 
of the first 150 day "long-run" field (LRcl) started in spring 2007. This field has been observed 
with BEST in order to identify variable stars. COROT is expected to obtain a precision of 0.66 
ppm in 5 days of obse rvation in the aste roseismic channel and 700 ppm in 1 hr of observation in 
the exoplanet channel (|Baglin et al.ll2002l ). This is an order of magnitude better than what can be 
obtained with BEST, because COPOT benefits from not having to observe through the atmosphere. 
Therefore, many of the transits that are expected to be detected with COPOT will not be visible 
in the BEST data, but despite this order-of-magnitude gap the BEST data will be highly valuable 
to the COPOT mission. 

In the exoplanet channel of the COROT satellite, stars have a large point-spread function 
(PSF) due to a prism installed in front of the detectors to obtain the colors of each star. The color 
information will help to distinguish stellar activity from exoplanet transits. In dense target fields, 
however, the resulting relatively large PSFs of the stellar spectra can overlap. Prior knowledge 
of variables overlapping the PSFs of interesting target stars will therefore be helpful to avoid 
misinterpretation. 



In addition to searching exoplanetary transits, additional science programs will be performed in 
the COPOT fields on selected targets. Some of these additional programs are interested in variable 
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stars, e.g., binaries. To be able to place integration windows in the satellite field on such targets, 
their position has to be known before observing the field with COROT. A search for variable stars 
in the COROT fields is therefore also helpful for the additional science program of the mission by 
detecting previously unknown variables in the field and confirming or rejecting the variability of 
the stars mentioned in variable star catalogs. 

A variability characterization like the one presented here of the target fields prior to COROT 
observations is therefore important for two reasons: to find interesting variable stars to monitor 
further by preselecting windows and to make an analysis of the impact of the variable stars on 
detecting extrasolar planet transits with COROT . 

The present paper is arranged as follows: Section 4.2 describes the observations. Section 4.3 
summarizes the data reduction procedure. Section 4.4 describes the data analysis, including search 
and detection of periodic variable stars. Section 4.5 contains the variable star catalog of the observed 
field, and Section 4.6 describes the analysis of the completeness of the data set. Concluding remarks 
are found in Section 4.7. 



2. Observation 

BEST is a small wide-angle telescope operated by the Institut fiir Planetenforschung of the 
Deutsches Zentrum fiir Luft- und Raumfahrt (DLR) and is located at the Observatoire de Haute- 
Provence (OHP), France. It consists of a Schmidt telescope with a 20 cm effective aperture, an 



f /2.7 focal ratio, and an AP-10 CCD camera, with 2048 x 2048 pixels of 14 fim size (iRauer et al 



20041 ). The result is a 3.1° x 3.1° field of view (FOV) with a pixel scale of 5'.'5. This is comparable 



to COROT, which has a pixel scale of 2'.'3 and a FOV for each of the four CCDs of 1.3° x 1.3°. 

We present the results from a monitoring program of the COROT field LRcl in the Galactic 
center direction centered at (a,5)=(19^00"^0a;0,+00°01'55'.'2) (J2000.0) obtained from 2005 June 6 
to September 7 . The data set consists of a total of 35 nights of observations spread out over 98 
nights in summer 2005. No data were obtained on nights with bad weather, around midsummer, or 
around full Moon. The observations were made without a filter in front of the CCD (the spectral 
response of the CCD corresponds to a wide R bandpass) in order to detect as many photons as 
possible and to provide a maximum sampling rate. An equal number of frames were obtained with 
40, 120, and 240 s integration times in order to observe stars in a large magnitude range and to 
analyze the effect of crowding in the frames. We only present data from the 240 s integration frames 
here, corresponding to 727 frames in total. Fig. 4.1 shows the 3.1° x 3.1° FOV observed. 
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Data Reduction 



The reduction of tiie data was performed as described by 



Based on a co mparative 

test of three different pho tometry al gorithms (S Extractor (jBertin &: Arnouta 119961 ) . MOMF 
(jKjeldsen k. FrandsenI Il992l ) and ISIS (jAlardI 120001 )) a data reduction pipehne was built (?) in 
order to reduce noise originating from the stellar background, extending the data pipeline of BEST 
data, which was optimized for less crowded fields observed at its previous location at the Thiiringer 
Landessternwarte, Germany (Rauer et al. 2004). We outline the different steps in the pipeline 
below. 



The frames a re calibrated from night to night by standard procedures as described in 



Rauer et al.ll2004l ). Hereafter, the reduction generally follows the method given in the ISIS package 



Alard 



20001 ). 



3.1. Image Subtraction 

Image subtraction is performed in order to remove both the sky and the stellar background. 
This is done by subtracting a reference frame from all the frames, which removes the sky background 
and all the constant stars. When image subtraction is performed in the right way it will produce 
uncrowded frames where only the variable stars are present. The noise will then be dominated by 
the photon noise from the target stars and from the stellar background. 

Before image subtraction can be made all the frames need to be transformed into the same 
(x, y)-coordinate grid. The frame with the best seeing from the middle of the observation run was 
selected as the template frame for this transformation. Then ISIS produced a reference frame by 
combining a number of best-seeing frames. The reference frame was then subtracted from all the 
other frames with the ISIS parameters listed in Table 4.1. This means that the reference frame is 
convolved with a convolution kernel in order to match all the stars in the frames. 

Having subtracted the reference frame from all the other frames, ISIS constructs a frame 
that is the mean of the absolute normalized deviations of all the subtracted frames. This frame 
only contains nonconstant stars, as all the constant stars have been removed by the subtraction. 
Unfortunately, this only works for nonconstant stars that have a large amplitude, as is discussed 
below. 

The ISIS package also includes a photometry routine that performs PSF-weighted photometry 
on the variable stars detected, thereby producing light curves of the differential fluxes of the stars. 
For different reasons, discussed below, we did not use this routine. We used the output of ISIS 
to distinguish between true periodic variable stars and stars that showed variability because of 
contamination from another nearby variable star. 
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3.2. Photometry 



There are three reasons why we do not use the photometry algorithm in ISIS. The first is that 
ISIS only produces light curves of differential fluxes; in order to convert these into magnitudes, 
photometry needs to also be performed on the referenc e frame in the same way as o n the subtracted 



frame s. The second reason, which has been noted by iMochejska et al.l (120021 ) and iHartman et al 



(|2004l ) ■ is that the detection of variables in ISIS works best for variables with an amplitude of 
at least 0.1 mag, since small-amplitude variables tend to get washed out in the overall noise in 
the subtracted frames. The last reason is that in order to do PSF-weighted photometry, the 
stars need to have a well-defined PSF; this is not always true if the stars are undersampled. We 
therefore preferred to perform simple unit-weight aperture photometry on the subtracted frames. 
Hartman et al.l (|2005l ) analyzed differences in performing PSF-weighted and unit-weighted aperture 
photometry on the subtracted frames and concluded that PSF-weighted photometry works best for 
the relatively faint stars, while unit-weighted aperture photometry works best for the relatively 
bright stars, but also that PSF-weighted photometry gives rise to a constant error term in the light 
curves. In addition, it is important that the transformation from differential fluxes to magnitudes 
is correct when evaluating the rms noise level obtained by different photometry algorithms and 
when comparing the noise level for ground-based transit surveys like BEST to that of space-based 
transit surveys like COROT. By using the same unit-weighted aperture photometry algorithm on 
the subtracted frames, as well as on the reference frame, we are confident that this transformation 
is correct. 



Hartman et al.l (|2004l ) note that if they use an aperture radius below 7 pixels for stars with 



a FWHM of the order of 3 pixels, the amplitude of the light curves is artificially reduced. As the 
FWHMs of the stars we analyze here are in the range from 1 to 3 pixels, depending on the seeing 
on different nights, we use an aperture radius of 7 pixels. We did try to reduce the aperture radius, 
which resulted not only in reduced amplitudes of the light curves but also reduced rms noise level 
in the light curves, as expected because of a decrease in the amount of background noise. In other 
words it might be possible to increase the precision by using smaller or nonconstant apertures, but 
then the amplitude and noise determination would not be as reliable. It is also noted here that 
although the PSFs in the nonsubtracted frames may be small, there may be large non-Gaussian 
residuals in the subtracted frames that require large apertures in order to ensure flux preservation. 
It is important to note here that the aperture radius of 7 pixels is only used on the subtracted 
frames where all the constant stars are removed. Therefore, using an aperture radius of 7 pixels 
does not mean that we are not able to resolve stars separated by less than 7 pixels (as long as they 
are not variables). Figure 4.2 shows the rms scatter in a single night (June 29) of all the light curves 
as a function of stellar luminosity, with the identified periodic variable stars marked. The figure 
shows that periodic variable stars are found roughly equally distributed over stellar luminosity. The 
obtained precision is comparable to what has been obtained by other small wide-angle telescopes 
([Hartman et aLlbooi IXane et aklbool boOSal lbll. 
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3.3. Astrometry and Calibration 



As noted bv lHartman et al.l (j2005l ) the light curves need to be cahbrated for systematic changes 
in the zero points. This is done by selecting a number of nonsaturated stars with low rms noise in 
the light curves and then calculating the mean change of the light curves. We selected 2182 stars 
for calculating these zero-point changes and found that the corrections from frame to frame were 
below 0.02 mag, except for a few bad frames which contained clouds or airplane tracks. 

Even after the zero-point correction, residuals of systematic trends in the light curves caused 
by, e.g., air mass, cloud cover, and temperature variations were still present in the data (these 
res iduals are sometime s referred to as "red noise"). We therefore applied the algorithm introduced 
by iTamuz et al.l (|2005l ) in order to correct for systematic effects. This clearly improved the rms 
noise level, in particular in the brightest stars. A detailed discussion of the photometric limitations 
of the BEST survey will be presented in a forthcoming paper. 

The astrometric calibration was performed by matching the 1500 brightest stars in the stellar 
catalog obtained from t he reference frame with the 2000 brightest stars (magnitude limit 11.4) of 
the USNO-A2.0 catalog (jMonet et al.lll998l ). The problem of matching the observed stars to stellar 
catalogs is not trivial. The main reason for this is the large image scale of BEST (5'.'5 pixel"^. 
This means that stars in a crowded part of the frame will have m ultiple matches in the stellar 
catalogs. The astrometric calibration was performed using MATCH (IValdes et al.lll995l ). MATCH 
finds a linear, quadratic, or cubic transformation between two stellar catalogs. We used a cubic 
transformation to calculate the transformation from the {x, y)-coordinates of the CCD to right 
ascension and declination. We note that a linear transformation does not work because of the large 
FOV. 

Out of the brightest 1500 stars in the stellar catalog obtained from the reference frame it was 
possible to match 902 with a matching radius of 2-10~^ degree. From these 902 stars we calculated 
an absolute magnitude difference between the observations and the USNO-A2.0 catalog R band. 
We note that we do not intend to perform a real absolute photometric calibration of the stars, since 
only relative magnitudes are needed to identify variables. 



4. Data Analysis 



4.1. Identification of Suspected Variable Stars 



We used the variability index J defined bv lStetsonI (|l996l ) to identify stars as suspected variable 
stars. The index J is defined as the normalized sum of the deviation of each data point (or pair of 
data points with small separations in time) in the light curve compared to the expected noise level 
for that star (the expected noise level is calculated based on the signal-to-noise ratio of the star). 
We considered data points to be pairs if they were separated by less than 0.03 day. Fig. 4.3 shows 
the index J as a function of magnitude. 
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We classified all 9112 stars with an index J greater than 2 as suspected variables, as all 
these stars show some kind of variability. Most of these stars are clearly variable stars, but either 
the period cannot be properly determined because of an imperfect duty cycle or the stars are 
nonperiodic or multiperiodic variable stars. In this paper we concentrate on the periodic variable 
stars. 



4.2. Identification of Periodic Variable Stars and Period Determination 

All the s tars classified as suspec t ed va riables were searched for periodicity with the method 



introduced by ISchwarzenberg-Czernvl (|l996l ). This method fits a set of periodic orthogonal polyno- 
mials to the observations and evaluates the quality of the fit with the use of an analysis-of- variance 
statistic. All the suspected variable stars were fitted with two polynomials, and all objects with a 
fit quality better than 0.9 were examined visually. 

In the visual examination we selected stars as periodic variable stars if they showed a clear 
periodicity between 0.286 and 120 days. We rejected light curves with periods close to 1 day, or 
harmonics of 1 day, because the periodicities in these light curves are likely to be caused by periodic 
changes in background level, temperature, air mass, etc. 

We have identified 92 periodic variable stars in the data set (Tables 4.2 and 4.3). It is seen 
in Fig. 4.1 that the periodic variable stars are detected all over the CCD, although the density 
depends on the stellar density in that part of the field. 



5. Variable Star Catalog 



The GCVS (iKholopov et al.l Il998l ) contains 60 variable stars in the field, 11 of which have 
periods less than 120 days. We have identified all 11 stars in the BEST data set. The phased light 
curves are shown in Fig. 4.4, and the corresponding parameters are given in Table 4.3. By looking 
at the light curves and comparing them with the index J we were able to confirm variability in six 
cases (Gl, G2, G8, G9, GIO, and Gil). None of the 11 stars were detected by the algorithm we had 
set up to search for periodicity. G8, G9, and Gil are RR Lyrae stars with such a characteristic light 
curve that they should have been detected by the algorithm. By going back in the analysis to see 
why these stars were not detected, we see that they were all excluded durin g the visual inspection 



There are three reasons for this exclusion. First, the method introduced by ISchwarzenberg-Czerny 



(|l996l ) provides a period that is equal to half the true period. This means that the phased light 
curves mainly show a falling trend. Second, the method introduced by Schwarzenberg-Czerny 
(1996) provides a period offset from the true period (the period of G8 is off by 5%, G9 by 0.2%, 
and Gil by 33%), thereby making the phased light curve look a lot more noisy. Third, G8 and 
Gil have periods very close to 1 day. Gl, G2, and GIO do not show strongly periodic signals, and 
detection was not expected. 
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We were not able to confirm variability in five cases (G3, G4, G5, G6, and G7). There are 
three possible reasons why we were not able to confirm variability in these three stars: (1) the 
identification as variable stars in the GCVS is wrong, (2) the stars are not variable anymore, and 
(3) the variability of the stars is below the noise level in the BEST data. In Table 4.3 it is seen 
that all three stars should have amplitudes larger than 1 mag, and in Fig. 4.4 it is seen that the 
observed variability in the light curves is well below 1 mag. This rules out the third possibility. 

We have divided the periodic variable stars identified by our search into two groups: periodic 
pulsating variables and eclipsing binaries. The two groups each contain 43 stars. For most of the 
stars it is easy to make this classification, but for some stars the classification may not be unique. 
The light curves of the periodic variable stars are shown in Figs 4.5 and 4.6, and the parameters 
are listed in Table 4.2. 

The detected periodic variable stars were cross-identified with the USNO-A2.0 catalog. For 
each star we checked that the cross-identification was unique. If there were multiple or no identi- 
fications we did not assign a USNO-A2.0 ID to the star. The USNO-A2.0 IDs are given in Table 
4.2. We also looked up all the detected periodic variable stars in the SIMBAD catalog but were 
only able to identify one of the periodic variable stars, i.e., V66 = V1225 Aql. 

Some of the stars show problems with blending. This means that the stars next to the variable 
stars also show variability. For each case where we had two or more stars showing the same kind of 
variability within the radius of the aperture (7 pixels) we have only identified one of the stars as a 
periodic variable. It is difficult to verify from the light curves alone which star is the true variable, 
because a faint star next to a bright variable star can show larger relative variability than the true 
bright variable star. In order to overcome this problem we identified all the variable stars that were 
affected by blending in the variability frame of our ISIS pipeline (the frame that was the mean of 
the absolute normalized deviations of all the subtracted frames), because this frame will only show 
truly variable stars. This approach only works well if the relative amplitudes in the variable stars 
are high. 

A few stars, such as V20, V31, and V50, seem to show some additional variability apart 
from the periodic signal. By inspecting the light curves we found that the additional variability 
originated from a few nights where the noise level was unusually high. It is therefore most unlikely 
that this additional variability is of stellar origin. 

Forty of the 92 identified periodic variable stars are in the LRcl FOV of one of the two COROT 
CCDs dedicated to the planet program. These are V3, V4, V6, V9, Vll, V12, V13, V15, V21, 
V23, V25, V26, V27, V32, V38, V39, V40, V41, V42, V43, V44, V45, V47, V49, V50, V53, V54, 
V58, V59, V60, V61, V65, V69, V78, V83, V84, V85, V86, G2, and G9. Some of the remaining 
identified periodic variable stars are expected to be in the FOV of some of the future pointings. 



We have identified all the periodic variable stars in the ExoDat catalog (IDeleuile et al.ll2006l ). 
The ExoDat catalog is part of the preparation of COROT and has been made in order to select 
which stars to observe with COROT, taking into account spectral type, luminosity class, variability. 
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and contamination (overlapping PSF on the focal plane). It contains, among other things, BVRI 
photometry of the field observed with the 2.5 m Isaac Newton Telescope at La Palma. 

We cross-identify the periodic variable stars from BEST with the ExoDat catalog using a 
matching radius of 1-10~^ deg. In this way we obtain 81 matches out of the 86 periodic variable 
stars. More identifications could have been obtained by increasing the matching radius, but that 
would have also increased the number of mismatches. 

Fig. 4.7 shows a color-magnitude diagram of the ExoDat catalog, with the discovered periodic 
variable stars marked. The stars in the V - I versus / diagram follow two branches: to the left 
the main-sequence stars and to the right the giant stars. For the faint stars there is a mismatch 
between stars on the main sequence and stars on the giant branch, since we do not know anything 
about the distance to the stars and therefore cannot provide their absolute luminosity. There is also 
a mismatch for nearby late main-sequence stars, which are bright with red colors and fall within or 
to the right of the giant branch. The diamonds in the figure represent eclipsing binaries, and the 
crosses represent periodic pulsating stars. A vague tendency for eclipsing binaries to fall mainly on 
the main sequence can be seen (15 out of 43; 35%), compared to the number of periodic pulsating 
stars that fall on the main sequence (7 out of 49; 14%). This is in agreement with what we expect, 
since most of the periodic pulsating stars we have found are Cepheids, Miras, or RR Lyrae stars, 
and none of these stars fall on the main-sequence branch in the V - I versus / diagram. 

Due to their short periodicity, the identified eclipsing binaries cannot be comprised of giant 
stars. Hence, from their location in Fig. 4.7, they range from solar-type stars (at the typical 
main-sequence field star distance of COROT Helds) to later but more nearby main-sequence stars. 



6. Data Completeness 



A way to evaluate the capability of ground-based surveys to detect variable stars is to compare 
the number of eclipsing binaries detected with a certain system with the expected number of 
eclipsing binaries in general. Evaluating the capability of ground-based surveys this way has the 
advantage that it includes not only the photometric precision of the data in the evaluation but also 
the phase coverage and the detection algorithms. 

Fig. 4. 8 shows the number of eclipsi ng binaries in the BEST data compared with the HIPPAR- 
COS data (ISoderhjelm &: Dischlen l2005l ) . A description of how the eclip sing binaries are detected 
in the HIPPARCOS data can be found in (jSoderhielm fc Dischleilbond ). 



The HIPPARCOS stars have been separated into two different absolute magnitude bins (My 
= -2.2 to 0.9) and (M^ = 0.3 to 4.0). This s eparation corresponds to a se paration into two different 
parallax bins ((vr > 4.5 and vr = 1.0 - 4.5) ISoderhjelm &: Dischleii l2005l ) . As we do not know the 
parallaxes for the stars observed with BEST, we do not know the absolute magnitudes. We have 
tried to compare the BEST data with the HIPPARCOS data from both bins with more or less the 
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same result, although the close HIPPARCOS sample seems to fit the BEST data slightly better 
than the distant sample. 

Fig. 4.8 shows that BEST detected more eclipsing binaries at low periods, which is caused 
by the higher sampling rate of BEST compared to HIPPARCOS. The fact that BEST detected 
more eclipsing binaries than HIPPARCOS at low periods indicates that HIPPARCOS may not 
be perfect at low periods. At higher periods (log P = 0.0 to 0.5 ) BEST only detected 20% - 



30% of the number of eclipsing binaries detected by HIPPARCOS. ISoderhjelm &: Dischleii (|2005l ) 
have calculated the relative number of eclipsing binaries as a function of period in three different 
amplitude bins (Am > 0.1, 0.4 and 0.7 mag). As almost all the eclipsing binaries we have detected 
have amplitudes between 0.1 and 0.4 mag, we only discuss the results from the first bin here, 
i.e., the bin with the largest eclipse probabilities. This statistical analysis shows that the BEST 
data set presented here has a completeness of 20% - 30% compared with the HIPPARCOS data 
set for periods longer than 1 day and is complete compared with the HIPPARCOS data set for 
short-period binaries. 

As there were six periodic variable stars in the GCVS that we did not detect with our periodic 
variable star detection algorithm, there is reason to believe that our algorithm is not perfect. This 
means that by improving the algorithm more periodic variable stars would be detected and the 
completeness of the data set would increase. In particular, it is seen that we detect a lower number 
of eclipsing binaries with periods around 1 day than what is expected. This might mean that we 
have been too conservative with excluding periodic variable stars that have a period close to 1 day. 



Conclusions 



BEST has observed the first long-run field (LRcl) of the COROT mission for 35 nights. We 
have found 92 periodic variable stars in the data: 43 eclipsing binaries and 49 periodic pulsating 
stars. 

Data reduction has been performed using image subtraction. Astrometric and photometric 
cal ibration were made by matching stars to the USNO-A2.0 catalog. We used the index J defined 
by IStetsonI (jl996l ) to ident ify suspected variable s t ars a nd searched these stars for periodicity with 
the method introduced bv ISchwarzenberg-Czernvl (|l996l ). In this way we have identified 92 periodic 
variable stars, compared to the already known 11 variable stars in the GCVS with periods less than 
120 days. 

The periodic variable stars have been identified in the ExoDat catalog, and the color-magnitude 
diagram of the stars in the field shows that the eclipsing binaries have a slight tendency to lie on 
the main sequence compared to the periodic pulsating stars. This refiects that most of the periodic 
pulsating stars detected in the data are Cepheids, Miras, or RR Lyrae stars. 



We also present a new tool to evaluate the completeness of a survey data set by comparing 
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the number of eclipsing binaries detected in the data set with statistical estimates of the relative 

number of eclipsing binaries. The analysis shows that this BEST data set has a completeness of 
20% - 30% compared with the HIPPARCOS data set for periods greater than 1 day and is complete 
relative to HIPPARCOS for shorter periods. 
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Table 1: Parameters used in ISIS 



sub_x: 


10 


Sub_y 


10 


deg_bg: 


3 


deg_spatial: 


2 


ngause: 


3 






sigma_gaussl: 


0.7 


deg_gaussl: 


6 


sigma_gauss2: 


2.0 


deg_gauss2: 


4 


sigma_gauss3: 


4.0 


deg_gauss3: 


3 


nstampsjx: 


10 


nstamps_y: 


10 


half_mesh_size: 


9 


half_stamp_size: 


15 
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Fig. 1. — The reference frame. The newly identified periodic variable stars are marked (V) as well 
as the stars from the GOVS (G). 
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Fig. 2. — The rms scatter of all the light curves {gray dots) and of the identified periodic variable 
stars (crosses). 
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Fig. 3. — The variability index J vs. R magnitude. The gray dots represents all the observed 
stars while the crosses represents the periodic variable stars that we have identified. The 9112 stars 
above the line with a index J larger then 2.0 were selected as suspected variable stars. 
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Table 2: Periodic variable stars in the field. Rmean is the mean magnitude obtained with BEST, 
Aj. is the amplitude of the periodic variability in the stars (corresponding to the deeper dip for 
eclipsing binaries) and the variability type is given as either periodic pulsating variables (PPV) or 
eclipsing binaries (EB). 



ID 


USNO 


a(J2000) 


5{J2000) 


P(days) 


^mean 


Ar 


Type 




VI 


U0826. 


.15631497 


19*'23"' 


50? 63 


-00° 34 


'53{'9 


1.908 


12, 


,21 


0, 


,06 


PPV 


V2 


U0826. 


.16062830 


19''30"' 


46? 56 


-00° 26 


'26"0 


11.932 


13. 


,40 


0, 


,16 


PPV 


V3 


U0826. 


.15906821 


19''28"' 


09? 70 


-00° 28 


'21('7 


0.487 


13. 


,70 


0, 


,16 


EB 


V4 


U0826. 


.15803107 




29? 85 


-00° 28 


'59"0 


2.113 


11. 


,64 


0, 


,04 


EB 


VB 


U0826. 


.15951762 


19*' 28" 


54? 77 


-00° 26 


'13" 9 


79.478 


12, 


,66 


0, 


,09 


PPV 


V6 


U0826. 


.15828621 


19*' 26™ 


54? 41 


-(10° 25 


'32" 8 


0.331 


13, 


,07 


0, 


,09 


EB 


V7 


U0825. 


.15479782 




36? 79 


-00° 27 


'18" 2 


1.967 


11 


,17 


0, 


.06 


PPV 


V8 


U0S25. 


.16122557 


19''31'" 


43? 22 


-00° 17 


'50" 6 


67.898 


10 


,30 





.12 


PPV 


V9 


U()S25. 


.15871361 




35? 18 


-00°17'49''5 


48.124 


12 


.50 





.07 


PPV 


VIO 


U0825. 


.16065804 


19''30™ 


49? 50 


-00° 12 


'59"4 


6.555 


11 


.26 





,02 


PPV 


Vll 


U0825. 


.15987353 


19'' 29" 


31? 05 


-00° 11 


'33" 6 


96.016 


10 


.76 





(16 


PPV 


V12 


U0826. 


.15757192 


19'' 25" 


46? 21 


-00° 11 


'33" 5 


90.120 


12 


.80 





15 


PPV 


via 






19'' 27" 


19? 88 


-00° 00 


'52;'6 


59.590 


11, 


.87 


0, 


.15 


PPV 


V14 


U0900. 


.15449170 




14?35 


00° 02' 


33^' 6 


102.780 


11, 


,37 


0, 


,10 


PPV 


VIB 


U0900. 


.15370565 


ig^'so™ 


15?97 


00° 05' 


37^' 6 


2.273 


12, 


,61 


0, 


,03 


EB 


V16 


U0900. 


.14944293 


19'>24" 


36? 78 


00° 01' 


16^' 6 


32.344 


11, 


,01 


0, 


,29 


PPV 


V17 


U0900. 


.15533956 


jghgjm 


17? 06 


00° 08' 


47" 7 


1.902 


12, 


,22 


0, 


,06 


EB 


VIS 


U0900. 


.14773005 


19^ 22" 


02?23 


00° 00' 


30^' 4 


1.135 


12, 


.22 


0, 


.13 


EB 


V19 


U0900. 


.14801088 


19''22" 


27? 37 


00° 06' 


59" 9 


1.870 


13, 


.16 


0, 


.06 


EB 


V20 


U0900. 


.14803632 


19''22" 


29? 64 


00° 09' 


57"4 


0.975 


13 


.16 


0, 


.18 


EB 


V21 


U0900. 


.14963651 


19'' 24" 


52? 18 


00° 20' 


50" 3 


85.208 


12 


.94 





.09 


PPV 


V22 






19''21" 


18? 72 


00° 17' 


50" 8 


93.828 


12 


.51 





,12 


PPV 


V23 


U()9()(). 


.15317860 


19'' 29" 


34? 44 


00° 32' 


54" 9 


15.250 


11 


.77 





.08 


PPV 


V24 


U0900. 


.14875071 


19'' 23" 


36?07 


00°31' 


39^' 4 


0.678 


11, 


.91 


0, 


.04 


EB 


V2B 


U0900. 


.15102073 


jghjgm 


47? 78 


00° 34' 47" 2 


78.260 


12, 


.62 


0, 


,07 


PPV 


V26 


U0900. 


.15105221 


jghjgm 


50? 33 


00° 34' 63" 4 


77.640 


12, 


.09 


0, 


,11 


PPV 


V27 


U0900. 


.15227975 


j^ghjgm 


26?95 


00° 43' 


19^'9 


116.850 


10, 


.64 


0, 


,12 


PPV 


V28 


U0900. 


.14801012 


19''22" 


27? 34 


00° 38' 


02" 4 


0.670 


12, 


.46 


0, 


,02 


EB 


V29 


U0900. 


.14785145 


19''22" 


13? 39 


00° 40' 


37''7 


43.136 


10, 


.71 


0, 


,04 


PPV 


V30 


U0900. 


.15459522 


19''31" 


21? 44 


00° 53' 


04^' 3 


66.864 


11, 


.70 


0, 


,07 


PPV 


vai 


U0900. 


.15338939 


jgh 2giii 


50?48 


00° 55' 


0l"6 


4.420 


12 


.28 


0, 


.04 


EB 


V32 


U0900. 


.15147997 


19''27" 


23? 94 


00° 53' 


07" 7 


0.454 


13 


.38 





.08 


EB 


V33 


U0900. 


.14857685 


19'' 23" 


20? 96 


01° 50' 


48" 6 


0.324 


14, 


64 





,19 


EB 


V34 


U0900. 


.14842306 


19'' 23" 


06? 29 


01° 06' 


49" 3 


0.549 


12 


.71 





(16 


EB 


V35 






19''31" 


48? 52 


01° 14' 


4l"7 


0.306 


13 


.72 





.13 


EB 


V36 


U()9()(). 


.15489834 


19'' 31" 


42?44 


01° 14' 


44" 2 


67.996 


12 


,38 





,07 


PPV 


V37 


U0900. 


.14842441 


19'' 23" 


06?42 


01°07'06^'8 


1.174 


12, 


,82 


0, 


,09 


EB 


V38 


U0900. 


.15096269 


^ghjgm 


42?91 


01°ll' 


48" 7 


0.452 


12, 


,20 


0, 


,07 


EB 


V39 


U0900. 


.15088039 


jgh2gm 


36? 27 


01° 13' 


IS^'S 


0.942 


13, 


,63 


0, 


,10 


PPV 


V40 


U0900. 


.15245627 


j^ghjgm 


37? 33 


01° 13' 


33^' 9 


1.783 


14, 


,66 


0, 


,10 


EB 


V41 


U0900. 


.14903885 


igh24" 


02?05 


01°2l' 


50^' 2 


39.304 


12, 


,18 


0, 


,20 


PPV 


V42 


U0900. 


.15033451 


19^ 25" 


51=33 


01°29' 


12" 5 


0.431 


14, 


,02 


0, 


,18 


EB 


V43 


U0900. 


.15172322 


19''27" 


43? 57 


01° 32' 


23"8 


63.682 


11, 


,62 


0, 


,03 


PPV 


V44 


U0900. 


.15110355 


jghj^m 


42? 76 


01°32' 


50" 7 


1.942 


11 


,37 


0, 


,03 


EB 


V45 


U()9()(). 


.15171348 


jghjgm 


54=39 


01°32' 


32" 1 


1.110 


13 


,23 


0, 


,10 


PPV 


V46 


U()9()(). 


.15245627 


19'' 28" 


40? 54 


01° 34' 


5()('4 


1.966 


11 


53 





,06 


EB 


V47 


U()9()(). 


.15190277 


19'' 27" 


58? 26 


01° 35' 


()4"l 


1.439 


11 


,61 





,07 


EB 


V48 


U()9()(). 


.15295393 


19'' 29" 


18? 23 


01° 37' 


29" 9 


55.496 


12 


99 





10 


PPV 


V49 


U()9()(). 


.149258()(i 


19'' 24" 


21? 19 


01° 33' 


31^7 


0.334 


12 


98 





(14 


EB 


V50 


U0900. 


.14877836 


19'' 23" 


38? 36 


01° 34' 


44" 


1.383 


12 


,35 


0, 


.06 


EB 


VBl 


U0900. 


.15318544 


^ghjgm 


34?91 


01°42' 


66'' 4 


62.763 


11, 


,17 


0, 


,06 


PPV 


VB2 


U0900. 


.14757657 


19''21" 


47? 43 


01°37' 


16;'3 


0.414 


13, 


,91 


0, 


,18 


EB 


VB3 


U0900. 


.14975844 


jghjgm 


02?35 


01°41' 


43;'0 


0.355 


11, 


,83 


0. 


,04 


EB 


V54 






^ghjgm 


59? 33 


01°4l' 


30^'l 


1.110 


13, 


,86 


0, 


,18 


PPV 


V55 


U0900. 


.15323691 


j^ghjgm 


38? 78 


01°46' 


35"4 


94.240 


11, 


,99 


0, 


,19 


PPV 


V56 


U0900. 


.14715405 


19''21" 


05?86 


01°46' 


4l''9 


89.746 


11, 


,80 


0, 


,06 


PPV 


V57 


U0900. 


.15337670 


jghjgm 


49?43 


01° 54' 


27"4 


38.362 


12 


,10 


0, 


,02 


PPV 


VB8 


U0900. 


.14845601 


19'' 23" 


09=52 


01°49' 


30" 9 


47.360 


10 


,07 


0, 


,10 


PPV 


VB9 






19'' 26" 


37? 45 


01°52' 


57^' 2 


0.286 


13, 


.11 


0, 


.08 


EB 


V60 


U0900. 


.15090051 


jgh2gm 


37?87 


01° 56' 


02" 9 


90.288 


12, 


,02 


0, 


.10 


PPV 
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VI 


P= 1.9080 






0-5 


1 ,0 

P=1.9665 




- 


0.5 

VI 


1.0 

P = 6.. 555.3 






0.5 

VI 3 


1.0 

P = 59.590 


- i.' 

i : ' 


1 ! ll 

y ^ 


0.5 

V21 


1.0 

P = 85.208 




0.5 

V25 


1.0 

P = 78.260 
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Fig. 4. — Light curves of 43 periodic pulsating stars in R. Magnitude is shown versus phase. The 
period is given in the upper right of each diagram is days. 
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0.5 


1-0 


13.03 


V45 




P=1 .1098 








13.16 


. 




• , ,J;'*v; • ■ 


13.29 








1 3.43 


_ 








0.5 


1.0 


13.45 


V54 




P=1 .1 lOJ 








13.72 
13.98 








1 4.25 










0.5 


1.0 


12.002 


V57 




P = 38.362 








12.069 


lili 

m 




'1 Pit: 


12.135 








12.202 








V62 


0.5 


1.0 

P = 74.072 


12.09 








12.26 


-• ■ 1 






12.43 


ll 




12.59 








V66 


0.5 


1.0 

P = 94.194 



0.5 1.0 



P = 55.496 V51 

1 1 .0 I 



I, 



P=94.2-10 Va6 





P=47.360 V60 







1,82 






' . i 




' 1 .95 
'2.09 


ft'' 

•f • '• 


•/ 






12.22 






0.5 
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0.5 
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Fig. 4. — Continued 
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Fig. 4. — Continued 



Table 2: Continued. 



ID 


USNO 


a(J2000) 


(5(J2000) 




P{days) 






An 


Type 


V61 


U()9()()_15089049 




37? 11 


01° 56 


17" 





2.045 


12 


02 


0.06 


EB 


V(i2 


U()9()()_1477()477 


19''22" 


02? 70 


01° 53 


01" 




74.072 


12 


34 


0.11 


PPV 


V63 


U0900.15513261 


19'' 32" 


00? 07 


02° 04 


19^' 


9 


51.242 


11 


60 


0.06 


PPV 


V64 


U0900.1B2864B7 




11?40 


02°02 


07^' 


5 


75.340 


11 


92 


0.06 


PPV 


V6B 


U0900.1B129182 


19''27'" 


08? 90 


02°04'04y 


2 


0.349 


12 


59 


0.05 


EB 


V66 


U0900_14931322 


19'' 24"" 


2B?96 


02°01 


51" 


4 


94.194 


11 


70 


0.08 


PPV 


V67 


U0900.14809956 




35? 19 


02° 02 


04" 


5 


1.362 


12 


87 


0.09 


EB 


V68 


U0900.15529595 




13? 30 


02° 12 


33" 


4 


73.852 


10 


69 


0.08 


PPV 


V69 


U()9()()_15083048 


19'' 26"' 


32? 09 


02° 07 


48" 


6 


76.042 


11 


40 


0.03 


PPV 


V70 


U0900_14819436 


19'' 22"' 


44? 20 


02° 05 


24" 


1 


93.784 


11 


02 


0.10 


PPV 


V71 


U0900.15313364 


19'' 29"' 


31?49 


02° 12 


54" 


6 


0.655 


13 


92 


0.05 


EB 


V72 


U0900.15131911 


19'' 27"" 


11?01 


02° 11 


09" 


9 


0.357 


13 


10 


0.14 


EB 


V73 


U0900.14936196 


19'' 24" 


30? 18 


02° 12 


53" 


3 


0.456 


13 


26 


0.18 


EB 


V74 


U0900.15290990 


19'' 29" 


14? 99 


02° 18 


20" 


5 


0.588 


13 


77 


0.10 


EB 


V7B 


U0900.14927833 


19'' 24"" 


22? 78 


02° IB 


02^' 


3 


1.028 


11 


62 


0.01 


EB 


V76 


U0900.1B07671B 


19''26" 


26? 80 


02° 18 


43^' 


4 


1.793 


10 


07 


0.02 


EB 


V77 


U0900.1B286206 


jghjgm 


11? 29 


02° 21 


65^' 


2 


0.392 


11 


99 


0.06 


EB 


V78 


U0826.1B906838 


jgh2gm 


09?62 


-00° 21 


'B8"2 


0.49B 


14 


21 


0.4 


EB 


V79 


U0900_14772878 


19''22"" 


02?07 


00°01 


07" 


5 


0.723 


13 


39 


0.2 


EB 


V80 


U0900.14716806 


jghj^m 


07? 48 


00° 33 


09" 


9 


1.188 


13 


26 


0.07 


PPV 


V81 


U0900_15357089 


19'' 30" 


05?' 04 


00° 46 


28" 


7 


0.307 


13 


40 


0.1 


EB 


V82 




19'' 23" 


28? 00 


00° 41 


29" 


1 


0.395 


13 


04 


0.05 


EB 


V83 


U0900_14966595 


19'' 24" 


54? 54 


01°15 


03" 


3 


5.872 


12 


20 


0.04 


PPV 


V84 


U0900.150334B1 


19''25" 


49? 73 


01°28 


43" 


2 


0.,549 


13 


95 


0.12 


EB 


V8B 




j^gh 25 m 


50? 12 


01°29 


43" 


5 


0.354 


13 


95 


0.2 


EB 


V86 


U0900.1B116292 


19'' 26"" 


5B?97 


01°31 


47^' 





2.B04 


13 


11 


0.6 


PPV 



Table 3: Periodic variable stars in the field from GCVS. 



ID 


GCVS 


c«(,72000) 


^(,72000) 


Typo 


P{days) 






flmoan 


Gl 


ES Aql 


19''32"21?50 


-00°ll'3l"8 


SR 


101.1 


13.2 


15.1 


10.80 


G2 


V0362 Aql 


19''24"55?67 


01°05'00"6 


RV 


68. 


13.1 


14.7 


14.09 


G3 


V0381 Aql 


19''31"13?68 


02°07'03''4 


RV 


54.8 


12.7 


15.4 


12.90 


G4 


V0919 Aql 


19'' 22" 55? 75 


01°02'39"7 


EB/KE 


0.797102 


13.0 


14.0 


13.98 


GB 


V0920 Aql 


19''31"30?76 


01°0B'53^'9 


RRAB 


0.4955033 


13.3 


14.3 


12.97 


G6 


V0291 Aql 


19'' 30"" 19? 93 


01°3B'47^'l 


RRAB 


0.482889 


13.6 


IB.O 


13.91 


G7 


V0922 Aql 


19''30'"28?36 


01°13'S0^'9 


RRAB 


0.4B8874 


16.3 


17.B 


12.99 


G8 


V0978 Aql 


19''31'"31?57 
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Fig. 5. — Light curves of 43 eclipsing binaries in R. See caption in Fig. HI 
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g. 5. — Continued 
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Fig. 5. — Continued 




Fig. 6. — BEST light curves of variable stars known from the General Catalogue of Variable Stars. 
See caption in Fig. [H 
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Fig. 7. — Color-magnitude diagram for stars in tlic COROT center filed. The gray dots are stars 
from the ExoDat catalog, the crosses are 49 pulsating periodic variable stars and the diamonds are 
43 eclipsing binaries found by BEST and identified in the ExoDat catalog. 
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Fig. 8. — Top: Relative number of detected eclipsing binaries as a function of period. The solid 
line shows relative number of eclipsing binaries i n the HIPPARCOS data - re lative to the total 
number of stars observed by HIPPARCOS (from ISoderhjelm &: Dischled (|2005l )) and the dashed 
line shows relative number of detected eclipsing binaries in the BEST observations - relative to the 
total number of stars observed in the field by BEST. Bottom: The relative difference between the 
relative number of detected eclipsing binaries in the BEST data set and in the HIPPARCOS data 
set - i.e. the relative difference between the two lines in th e upper plot. The HIPPARC OS data 
were obtained as the mean value of Fig. 14 and Fig. 4.15 in ISoderhjelm &: Dischleii (|2005l ). 



